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Abstract 

The Photogrammetric Appendage Structural 
Dynamics Experiment was designed, developed, and 
flown to demonstrate and prove measurement of the 
structural vibration response of a Russian Space Station 
Mir solar array using photogrammetric methods. The 
experiment flew on the STS-74 Space Shuttle mission 
to Mir in November 1995 and obtained video imagery 
of solar array structural response to various excitation 
events. The video imagery has been digitized and 
triangulated to obtain response time history data at 
discrete points on the solar array. This data has been 
further processed using the Eigensystem Realization 
Algorithm modal identification technique to determine 
the natural vibration frequencies, damping, and mode 
shapes of the solar array. The results demonstrate that 
photogrammetric measurement of articulating, non- 
optically targeted, flexible solar arrays and appendages is 
a viable, low-cost measurement option for the 
International Space Station. 

Introduction 

On-orbit structural response measurements of the 
International Space Station (ISS) are required for loads 
validation and verification of structural mathematical 
models. Typically, structural response time-history data 
would be collected with on-board sensors such as strain 
gages or accelerometers, signal conditioned, and 
transmitted (downlinked) to the ground for further 
processing. From the time history data, structural 
dynamic characteristics (structural mode frequencies, 
dampings, and mode shapes) would be determined using 
structural identification algorithms such as the 
Eigensystem Realization Algorithm (ERA) [1]. 
Accelerometer and strain gage measurements of the ISS 
primary truss are currently planned for loads validation 
and model verification; however, accelerometer or strain 
gage measurements of solar array response are not being 
considered because of the impact on cost and resources. 

Photogrammetric measurement of solar array 
structural response is a low cost alternative to dedicated 
accelerometer measurements, especially when 
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measurements are required for articulating or rotating 
spacecraft components such as solar arrays or thermal 
radiators [2,3]. Elimination of accelerometers and their 
associated power and data wiring, signal conditioning 
and digital conversion electronics, etc., can greatly 
simplify the spacecraft electrical design and integration 
and reduce spacecraft cost. 

The Photogrammetric Appendage Structural 
Dynamics Experiment (PASDE) was developed to 
mitigate the technical risk and cost associated with on- 
orbit measurement of solar array and other flexible 
appendage structural responses for the ISS program. 

The experiment objectives were to demonstrate 
photogrammetric structural response measurement 
methods for solar arrays using video cameras without 
the use of optical targets, to provide engineering data on 
solar array designs similar to those expected to be used 
on the ISS, and to verify that routine on-orbit spacecraft 
operational events provide sufficient excitation for 
structural response testing. 

The NASA Langley Research Center (LaRC) and 
the NASA Goddard Space Flight Center (GSFC) were 
funded to conduct the PASDE experiment. PASDE 
flew in space on the Space Shuttle STS-74 mission in 
November 1995 to the Russian Mir space station as 
part of the ISS Phase I program. For the STS-74 
mission, PASDE was manifested with the Arizona 
State University Glo-4 experiment as the Hitchhiker 
Glo-4/PASDE or GPP payload [4]. 

Design and Development 

The objective of the PASDE experiment was to 
measure the response of the Mir Kvant-II module lower 
solar array shown in Figure 1 . During the docked 
portions of the mission, this array was outboard of the 
Shuttle on the port side, slightly ahead of the wing. 
Based on the length of the Kvant-II solar array and 
possible viewing geometries, a PASDE design using 
six video cameras and recorders in three Hitchhiker 
canisters located in the Shuttle payload bay was 
selected. In this design, one camera in each canister is 
aimed and focused at the root of the array, and the other 
at the tip to provide full video coverage and optimize to 
the extent possible the triangulation geometry. Each 
video recorders could record 1 20 minutes of video data. 


Three sets of identical flight hardware were 
fabricated, assembled, and tested at the NASA Langley 
Research Center beginning in October 1 994. Each 
PASDE flight unit was integrated with a Hitchhiker 
canister at GSFC and the integrated Hitchhiker hardware 
was delivered to the NASA Kennedy Space Center 
(KSC) in July 1995. At KSC a pre-flight 
photogrammetric calibration of the hardware and ground 
simulation of the experiment was performed using a 
simulated solar array assembly. Results from this pre- 
flight ground experiment are documented in [5]. 

The three Hitchhiker canisters containing the 
PASDE flight hardware were mounted aft of the 
expected Kvant-H array position in Space Shuttle 
Atlantis payload Bay 6 (port side). Bay 7 (starboard 
side), and Bay 13 (port side) following the pre-flight 
calibration and testing. Additional information on 
PASDE design requirements and hardware are described 
in [6]. 

Mission Operations 

STS-74 launched at 7:30:43 am EST on November 
11, 1995. PASDE mission operations were conducted 
during the STS-74 mission from the Hitchhiker Payload 
Operations Control Center (POCC) at NASA GSFC. 
The PASDE hardware was commanded from the POCC 
by Hitchhiker personnel under the direction of the 
PASDE operations team. Table 1 summarizes actual 
PASDE operations during the mission. Events are 
listed with reference to Mission-Elapsed-Time (MET) in 
day: hours: minutes format, amount of data recorded in 
minutes: seconds format, and a short description. 
Additional details regarding the on-orbit operations of 
PASDE can be found in [6], 

Post-Flight Calibration and Flight Data 
Retrieval 

The PASDE hardware underwent post-flight 
photogrammetric calibration following removal from 
the Shuttle using the same procedures and simulated 
solar array as were used in the pre-flight calibration. 

The post-flight calibration data was recorded directly on 
the flight data tapes and used for subsequent 
photogrammetric data analysis processing. Upon return 
of the hardware to LaRC, the PASDE data tapes were 
removed from the flight video recorders and duplicated. 
Five of the six camera/recorder combinations performed 
as designed, providing excellent quality video imagery 
of solar array motions. The Bay 7 Tip camera exhibited 
intermittant anomalous behavior, resulting in some data 
loss (estimated at 50%) for that camera/recorder 
combination. 


Data Processing Methods 

The process of extracting the structural dynamic 
characteristics of mode frequencies, dampings, and mode 
shapes from the recorded video images is described 
below. The basic outline of the process is a) 
digitization (conversion) of the analog video data to 
digital form, b) identification and tracking of identifiable 
features in the digitized images for a single camera, in 
camera focal plane x-y pixel coordinates using line- 
tracking or scene correlation algorithms, c) triangulation 
of results for points coincident in 2 or more digital 
image sequences to physical (3 -axis) coordinates, and d) 
application of the Eigensystem Realization Algorithm 
to the time history data in physical coordinates to 
determine modal characteristics. 

Digitization - The video flight data was digitized 
using the hardware arrangement illustrated in Figure 2. 
Table 2 summarizes the data which has been digitized to 
date. The real-time video data rate was 30 frames per 
second. However, since the highest mode of interest is 
on the order of 0.5 Hz, a sampling rate of 1 0 fps was 
sufficient to characterize the motion of the array. 
Therefore, only every third frame of video data was 
saved for further processing. Figure 3 shows digitized 
images from five of the six PASDE cameras for the 
PRCS #3, Firing #2 data case (Note: usable data was 
not obtained from the Bay 7 Tip camera during this test 
due to the anomalous behavior previously mentioned). 

Line Tracking - In the absence of optical targets, 
naturally occurring lines and features in the digitized 
images may be tracked with respect to their x and y 
pixel locations using the method of line tracking 
described in [7]. As shown in Figure 4, two lines 
delineating a target of interest may be selected by a 2 
pairs of points. From the calculated slope and intercept 
of each line, their intersection point may be computed. 
This is the initial target position for the tracking 
algorithm. For subsequent images in a specified 
sequence, the tracking algorithm scans along each of the 
lines to find the peak pixel location (maximum 
intensity value) within a specified pixel width. The 
centroid of each scan is found, a straight line fit of the 
centroid values along each line is then made, and a new 
intersection point is computed. The x and y pixel time 
histories of the motion of an identified target are thereby 
generated for the sequence of images. 

Cross-Correlation Analysis - Correlation of images 
allows matching of particular features in two or more 
separate views [8]. The cross-correlation function 
measures the degree of similarity of two images at 
various amounts of shifting of one image relative to the 
other. Maximum correlation occurs at the x-y shift 
values where the two images are most similar. By 
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finding the point of maximum correlation, the change 
of position of a “target” (i.e., any identifiable feature) 
from one image to the other is determined- Applying 
this fact to a sequence of images generates a time 
history of the motion. 

Referring to Figure 5, the cross-correlation function 
of subimages A and B is calculated by shifting 
subimage A over the range of ± M pixels in the x 
direction and ± N pixels in the y direction from its 
center position. Subimage A (extracted from the overall 
measured image) is sized to contain characteristic 
contrast patterns in both the x and y directions centered 
on the “target” of interest. The shift values - M < m < 
M and - N < n < N must be somewhat larger than the 
distance that the target has moved from subimage A to 
subimage B. In the PASDE data analysis, a typical size 
of subimage A is 41 rows x 41 columns (from a total 
image of 480 rows x 512 columns), and a typical size 
of M and N is 4. 

In practice, cross-correlation calculations occur in 
one of two ways: 1) correlation of a fixed video frame 
with a sequence of succeeding frames, or 2) correlation 
of adjacent pairs of frames with summation of the total 
amount of shift. The first approach is computationally 
simpler but does not account for rotation or scale 
change of the image versus time. It is applicable for 
most vibration problems where motion oscillates about 
a static position of the structure. The second approach 
allows general tracking of motion and is applicable for 
problems with large changes of orientation such as solar 
array slewing. 

Figure 6 shows a typical cross-correlation function. 
Note that the basic correlation function, has a spatial 
resolution of 1 pixel. However the correlation function 
itself inherently has much higher resolution if other 
pixels in the neighborhood of the pixel of maximum 
intensity are considered. It is important in precision 
photogrammetry to interpolate the correlation function 
to subpixel resolution in order to locate as precisely as 
possible the position of maximum correlation [9]. The 
‘x’ in Figure 6 shows an interpolated position of 
maximum correlation for the typical case. The contour 
lines in Figure 6 show the general shape of the 
correlation surface but are not used directly in the 
interpolation calculations. 

Comparison of Target Tracking Methods - Both the 
cross-correlation and the line tracking method were 
applied to a target on the panel of the Kvant-II array as 
viewed from the Bay 13 Tip camera. Pixel time 
histories generated by each method are shown in Figure 
7. In this case, and in similar cases where there is a 
high amount of contrast in the image, the results from 
the two methods are very comparable. The slight 


difference in mean y pixel position is attributed to the 
fact that the initial targets were slightly mismatched. 
The frequency and trends in the time histories are in 
good agreement. 

Triangulation - The pixel time histories generated 
by the tracking procedures for identified targets in two 
or more cameras are triangulated to produce three-axis 
displacement of the target in the OSR coordinate 
system. The pixels are first mapped to inches in the 
image space coordinate system. These image space 
coordinates are then related to the object space 
coordinate system (OSR in this case) via the 
photogrammetric collinearity equations [3]: 
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where x and y are the coordinates of the image point, x Q 
and y G are the coordinates of the principal point, c is the 
principal distance of the lens, X, Y, and Z are the 
coordinates of the points in the object space coordinate 
system, X c , Y c , and Z c are the coordinates of the 
perspective center, k is a scale factor (= 1.0 here), and 
M is a 3x3 orthogonal rotation matrix the elements m^ 
of which are the cosines of the angles between image 
and object space coordinate systems. 

Determination of the (exterior and interior 
orientation) parameters of the collinearity equations was 
from the pre- and post-flight ground calibration data 
using the method of least squares. Using the computed 
camera parameters, triangulation of the time histories of 
targeted points from various camera views was 
accomplished. Defining 
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Equation ( 1 ) was rewritten for camera i as 


B,d=B,<i; 

where 


( 2 ) 

( 3 ) 


r(wii+*m 3 i) 

|_( m 2l +/"»3l) 


d = 


X 

Y , 
Z 


( m n + x'm i2 ) 

(m 22 + y'm 32 ) 


(m n + x'm 33 ) 
(m 23 + y'm 33 ) 


and 


3 




For a target viewed from l cameras, the vector of object 
space coordinates was computed by: 


d = N ‘t (4) 

where 


i i 

N = andt=£B'B,.d; 

1=1 J = 1 

Modal Identification - The objective of modal 
identification is to decompose measured free -decay time 
histories into their constituent structural modal 
components. Each mode can be characterized by three 
parameters: natural frequency, damping, and mode 
shape. The natural frequency and damping are scalar 
quantities, whereas the mode shape is a complex-valued 
vector of n components, where n is the number of 
discrete measurement points. The Eigensystem 
Realization Algorithm (ERA) [ 1 ] was used for PASDE 
modal identification. 

ERA analysis assumes that each of the 
measurements is the free-decay vibration response at 
some point on a linear, multiple-degree-of-freedom 
system. Using digitized samples of time histories 
obtained as described in sections above, x(j),j= 1,2,3,... 
ERA calculates the unknown parameters in the 
following characteristic equation: 

M 

x i (O = X A ir e ° r ' S ‘ n (®r* + <t>,r ) (5) 

r= 1 

where the unknown parameters are: M y the number of 
excited modes, A (r , the initial amplitude of mode r at 
measurement location /, <p in the initial phase angle of 
mode r at measurement location i , a) r , the frequency of 
mode r, and o r , the damping of mode r. 

With experimental data, the number of modes M is 
often ambiguous due to “real world” effects such as 
nonlinearities, noise, distortion, etc. Such ambiguity 
can cause significant fluctuations in some identification 
results. With experimental data the characteristic 
equation given above is not satisfied exactly for any set 
of parameters, and the best approximation (in a linear, 
least-squares sense) is used. When data for multiple 
measurements (and/or tests) are analyzed 
simultaneously, a single value of natural frequency and 
damping is obtained by ERA for each mode. Because 
natural frequencies and damping theoretically are global 
quantities (i.e., identical at all measurement locations 
and in all tests of the same structure), simultaneous 


analysis of multiple measurements often provides better 
ERA results than analysis of individual measurements. 

Results 

PRCS Test #3. Firing #2 - Structural dynamic 
characterization (determination of structural mode 
frequencies, dampings, and mode shapes) has been 
completed for the PRCS Test #3, Firing #2 (see Table 
1). Figure 3 shows 12 locations on the Kvant-II solar 
array selected for response time history calculations. 

The locations are equally spaced in 2 columns on either 
side of the centerline of the solar array. Location Nos. 

1 through 8 appear in all 3 sets of root camera video 
recordings. Location Nos. 9 through 1 2, however, 
appear in only 2 of the 3 sets of tip camera video 
recordings due to the data loss with the Bay 7 Tip 
camera previously noted. The motion at each location 
was determined using correlation analysis method, also 
as described previously. 

Figure 8 shows selected time history data obtained 
from the correlation tracking method from 4 of the 
twelve points identified in Figure 3 (points 4, 8, 10, 
and 12). The time histories are the final results of the 
photogrammetric analysis (triangulation) presented in 3- 
dimensional Orbiter Structural Reference (OSR) frame 
coordinates. The data are arranged in 4 rows by 3 
columns corresponding to the 4 measurement locations 
and the x, y, and z positions in the OSR. These data 
are of high quality considering the small amount of 
motion occurring in this test (a maximum of 1 inch 
peak-to-peak). The largest motions occur in the x 
directions and these data appear as classical free-decay 
vibrations over the majority of the time interval. There 
are some unexplained fluctuations and apparent drifting 
in the data. The most likely source of these deviations 
from classical free-decay vibration is friction (stiction) 
and/or freeplay in the mechanisms of the Shuttle-Mir 
mated configuration (including the Shuttle-to-Mir 
docking mechanism, the rotary drive of the solar array, 
and the hinges and motors used to unfold the solar array 
from its stowed position). The time histories for 
locations 4 and 8 are 1 40 seconds in length whereas 
those for Locations 10 and 12 are only 60 seconds in 
length. 

Five dominant modes were identified using ERA 
from the twelve points shown in Figure 3. Based on 
the computed mode shapes, shown in Figure 9, two of 
the identified modes were system modes of the mated 
Shuttle-Mir configuration, whereas the other three 
identified modes were flexible modes of the solar array. 
The two system modes have the following common 
characteristics: a) the largest motion appears at the root 
of the solar array rather than the tip, b) the direction of 
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motion is predominantly in the x direction (OSR) rather 
than in the solar array out-of-plane or in-plane directions 
(the solar array out-of-plane direction is rotated 25 
degrees from the OSR x axis), and c) the damping level 
is relatively high (23.7% for system mode 1 and 3.4% 
for system mode 2). Common characteristics of the 3 
solar array modes are: a) the largest motion appears at 
the tip of the solar array, b) the direction of motion is 
predominantly in the solar array out-of-plane (2 modes) 
or in-plane direction (1 mode), and c) the damping level 
is generally lower (2.4% for the first out-of-plane mode, 
1.7% for second out-of-plane mode, and 3.7 for the in- 
plane mode). 

Because basic ERA analysis is a linear modal 
identification technique, the results shown provide a 
good linear representation of the solar array structural 
dynamics. However, significant nonlinearities were 
also observed in this set of PASDE data. The nonlinear 
effects were characterized using a sliding time- window 
analysis (i.e., a piecewise linear analysis) within ERA. 
For a linear system, natural frequencies and damping 
factors are constant as a function of the time shift 
within the sliding time-window analysis. Various types 
of structural nonlinearities, on the other hand, produce 
various changes in frequency and damping as a function 
of time shift. 

The results of the sliding time-window analysis 
show that the frequency and damping are relatively 
constant for the first two solar array modes, while large 
fluctuations occur in the results for the first system 
mode. The first system mode probably involves the 
massive Mir structure flexing relative to the massive 
Shuttle structure through the flexibility (and damping) 
of the docking interface. This interface contains several 
mechanisms which account for the large observed 
damping and nonlinearity. Changes with time shift 
also occur for the other two solar array modes, 
indicative of appreciable nonlinearity, primarily in 
damping (perhaps energy exchange between modes). 

Table 3 further summarizes the modal identification 
results. The range of frequency and damping values are 
the maximum and minimum values observed in the 
sliding time-window analysis. The final 2 columns in 
Table 3 show the physical displacement amplitude of 
each mode. Both the maximum initial zero-peak 
amplitude and the root-mean- square (rms) amplitude in 
the data-analysis interval are reported. 

Night/Dav Terminator - The orbital night/day and 
day/night terminator crossings of the docked 
Shuttle/Mir spacecraft were of interest in PASDE 
because of potential thermally-induced structural 
responses of the solar array. Time history response 
analysis of the night/day terminator test (see Table 1) 


has been completed. Three-axis triangulated data for a 
point near the tip of the array (approximately location 
12 in Figure 3) is shown in Figure 10. Over a period 
of about 3 minutes following the entry of the combined 
Shuttle/Mir spacecraft into sunlight from darkness, the 
tip of the array moves quasi-statically about 0.6 inches 
in the x-axis direction combined with an approximately 
± 0.05 inch, 0.16 Hz oscillation. Additional analysis 
of data for points closer to the root of the array is needed 
to determine if the measured quasi -static response is 
limited to the array only, or is the response of the 
combined Shuttle/Mir spacecraft to thermal loading. 

Solar Array Slew - The solar array sun track 
rotations of the Kvant-II solar array were also considered 
to be a potential source of structural excitation because, 
unlike many U.S. systems, the Mir arrays track the sun 
in approximately 22 degree increments at relatively high 
rotational rates. Time history analysis of data for this 
test (see Table 1) is shown in Figure 1 1 (x and y axis 
results only). The triangulated results are for a point 
near the root of the solar array (approximately location 
2 in Figure 3). A solar array axial rotation angle (slew 
angle) has been computed from the data and is very 
close to the expected 22 degrees of rotation. Visible in 
the data during the slew is an 0.16 - 0.17 Hz oscillation 
corresponding to an excited structural mode. 

Concluding Remarks 

The Photogrammetric Appendage Structural 
Dynamics Experiment was flown on the Space Shuttle 
STS-74 Mission to the Russian Space Station Mir. A 
total of 1 1 3 minutes of video imagery of the Mir 
Kvant-II solar responding to various excitations was 
recorded. The data has been digitized, converted to three- 
axis time histories for various discrete points on the 
solar array, and analyzed using the Eigensystem 
Realization Algorithm to determine structural mode 
characteristics. Five structural modes were identified, 
two combined Shuttle/Mir modes, and three Kvant-II 
solar array modes. In addition, the response of the array 
to thermal loadings from night-to-day orbital transitions 
and solar array sun track slewing has been determined. 
The experiment demonstrated that low-cost, passive 
photogrammetric measurement of flexible solar array 
and appendage structural response is feasible for use on 
the International Space Station. 
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Table 1 - STS - 74 Mission Events 



EB9 

Description 

0:00:00 

00:00 

STS-74 launch at 7:30:43am EST, 11/12/95 

0:03:25 

00:00 

PASDE activation by Shuttle crew 

1:04:45 

11:35 

Camera alignment verification and contingency data collection 

2:17:52 

10:47 

Shuttle docking with Mir 

2:19:22 

18:34 

PRCS Test #1 

3:16:54 

09:30 j 

PRCS Test #2 

4:21:28 

02:00 

Day-to-night terminator (Bay 6 canister only) 

4:21:54 

05:14 

Night-to-day terminator 

5:00:35 

10:30 

Shuttle/Mir attitude maneuver using Shuttle primary jets 

5:01:43 

21:50 

PRCS Test #3 

5:06:15 

01:20 

Solar array sun track slew maneuver #1 

5:06:24 

08:35 

Solar array sun track slew maneuver #2 

- 

13:31 

Other data tape recording (recorder/commanding tests, etc.) 


113:26 

Total on-orbit record time (115 minutes available) 
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Table 2 - Digitized PASDE Flight Data Summ 


Event 


Night/Day Terminator 
PRCS #3 

F iring #2 

PRCS #3 

Firing #3 

Solar Array Slew 


MET start 


004/21:56:30.02 

005/01:54:00.02 

005/01:57:30.02 

005/06:34:50.02 


MET sto 


004/21:59:30.02 

005/01:56:20.29 

005/01:58:30.02 

005/06:25:30.02 





Nominal 

Nominal 

Max. Initial 

RMS 

Frequency 
Variation in 

Damping 
Variation in 

Mode 

Description 

i~~ i 

Frequency, Hz 

Damping 
Factor, % 

0-Pk 

Amplitude, in 

Amplitude 
in Analysis 
Window, in 

Sliding 

Window 

Analysis, 

Hz 

Sliding 

Window 

Analysis, 

% 



Figure 1 - STS-74 Shuttle/Mir configuration 
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Figure 2 - Digitization Layout 
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PRCS_3. Camero 5: ROOT 1 3-5. 1 .54.0.2.mat 
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Figure 3 - Locations on Kvant-2 Solar Array Centerbox Tracked With Correlation Analysis, 

Test PRCS #3, Firing #2 (+ Pitch) 
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Figure 4 - (Line Track Method) 
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Figure 5 - Correlation of Subimages A and B 


Figure 6 - Typical Cross-Correlation Function 
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Figure 7 - Comparison of Line Track (Solid) & 
Correlation Methods (Dash) 









Figure 8 - Displacements in Shuttle Coordinates at Locations 4, 8, 10, and 12 on Kvant-2 Solar Array Centerbox, 

Test PRCS #3, Firing #2 (+ Pitch) 
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Figure 9 - Identified Mode Shapes, Test PRCS #3, Firing #2 (+ Pitch) 
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Figure 10 - Displacements in Shuttle Coordinates Near Solar Array Tip, Night/Day Terminator Test 



Figure 1 1 - Displacements in Shuttle Coordinates Near Solar Array Tip, Solar Array Slew Test 
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